Negative differential resistance in nanoscale transport in the Coulomb blockade 

regime 
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Abstract 

Motivated by recent experiments, we have studied transport behavior of coupled quantum dot 
systems in the Coulomb blockade regime using the master (rate) equation approach. We explore how 
electron-electron interactions in a donor- accept or system, resembling weakly coupled quantum dots 
with varying charging energy, can modify the system's response to an external bias, taking it from 
normal Coulomb blockade behavior to negative differential resistance (NDR) in the curent-voltage 
characteristics. 



The switching and negative differential resistance (NDR) behavior of nanoscale systems has gained a lot of interest 
in the last decade, owing to the potential applications in single molecule electronics and has been observed in a variety 
of experimental systems, especially in the widely studied Tour molecules [J, Q|. There have been many theoretical 
studies to understand this phenomenon within mainly through the one-electron picture 3j, \% [Si Ifl 13] • There have 
also been a number of theoretical studies on donor-acceptor double quantum dot systems, where strongrectification 
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1J or 



has been observed [12J, and others which showed NDR with variation in the dot-electrode coupling [8|, 
due to a detuning of the dot levels [13| . Another recent study has attempted to establish the conditions obeyed by 
the parameters involved, to find such a collapse in the current magnitude . 



quantum dots also showed an NDR feature 



1G 



15 1 . Some recent experiments on double 



17] and has rekindled interest in the phenomenon, occuring in the low 



temperature, weak-coupling limit. Theoretical studies of NDR in this single electron charging limit, is now gaining 



prominence and attracting a lot of research 9 



10 



11| . This regime, where mean- field descriptions usually fail, is 



one where electron charging energies are very high as compared to the broadenings due to average coupling, and 



are particularly important for small molecules which behave more like a quantum dot than a wire |18l. 
mean field methods combined with standard non-equilibrium Greens function (NEGF) 



22, 



19(. Since 



231 treatment of 



transport is perturbative in the interaction parameter, it cannot capture the transitions between the spectrum of 
neutral and excited states, which can lead to a variety of interesting features in the current-voltage characteristics. 



* Presently, at Institute for Materials Science and Max Bergmann Center of Biomaterials, 
Dresden University of Technology, D-01062 Dresden, Germany 




FIG. 1: A schematic representation of a two-dot system consisting of a donor and an acceptor coupled to two electrodes. 
The formalism that has now come to be used widely to capture molecular transport in the Coulomb blockade regime 



is the master or rate equation method 



24|. 



In this article, we use the above formalism to study a two-dot system consisting of a donor and an acceptor (see 
the schematic given in figure 1) in the Coulomb blockade regime. Taking cue from our previous mean-field transport 
studies on two-level systems which showed interesting non- linear behavior in their current- voltage characteristics 6( , 
here we explore the role of strong correlations in affecting their transport behavior. This study becomes interesting, 
especially in the context of the difference in their low-lying excitations, which would play a very important role in 
their low-bias current- volatage characteristics. The rate equation formalism describes transport through a correlated 
system with many-body eigenstates. The presence of Coulomb interactions results in occupation probabilities of each 
many body state that cannot be factorized as the product of the occupation probabilities of each single electron level. 
Hence, in this case, the full rate-equation problem, where the occupation probability of each many-body state is 
treated as an independent variable is solved, neglecting off-diagonal coherences. In this method, the transition rate, 
S s /^ s , from the many-body state s' to s, differing by one electron, is calculated up to linear order in T (which is the 
bare electron tunneling rate between the system and the electrodes), using Fermi's golden rule as 



^,t s = rf L (E s -Ei s )J2\<s\clW>\ 2 

a 

= rfn(E s - E' S )J2\ < *\CLW > I 2 (1) 

a 

with a corresponding equation for and S^ s , obtained by replacing Jl,r(E s — E' s ) by (1 — fL.n(E s — E' s )). 

Here, +/— correspond to the creation/annihilation of an electron inside the dot due to electron movement from/to 
left (L) or right (R) electrode. C\ a and C\j a are the creation operators for electrons of spin, a at the first and Nth 
sites respectively. We have also assumed that the creation and annihilation happen only at the terminal sites. The 
total transition rate is then obtained as, S s ^ s ' = S^ s , + S^ s , + + £^ ,. The non-equilibrium probability 

P„ of occurrence of each many-body state s is obtained by solving the set of independent rate equations defined 



3 

by P s — ^ S ,(£ S /_> S P S / — Es—^'Pg) through the stationarity condition P s = at steady state. This results in a 
homogeneous set of equations of the size of the many-body space. Taking advantage of the normalization condition 
^2 S P S — 1, we obtain linear equations, which can be solved using well-known linear algebraic methods. The steady 
state probabilities are then used to obtain the terminal current as, 



tX-''-'' - ( 2 ) 



-S..S' 



where a = L/R. Using the above prescribed method, we study a two site system described by the Hamiltonian, 

2 2 

H = (ej - eW g )a\ai + -t(a\ a a 2(7 + h.c.) + U n^mi + Vu(ni - n)(n 2 - n) (3) 
»=i o-=T,l *=i 

where t is the hopping strength between the sites with same spin (tr), e\ t 2 are the on-site energies, U is the Hubbard 
interaction between electrons at the same site, V12 is the nearest- neig hbor Coulomb interaction and W g is the external 



gate bias. The average charge (n) is assumed to be unity here 26|. n actually gives a constant shift to the energy 
levels with fixed number of electrons. For two sites with two electrons, the energy levels are negatively shifted by 
V12 amount. Note that, there exists two quantum phases in this model with variation of interaction parameters. For 
the half filled ground state, with zero onsite energies, while U > V12/2 represents a spin, density wave (SDW) phase, 



27 , |28|1 . However, in our case with 



U < V12/2 corresponds to charge density (CDW) phase in thermodynamic limit 
two sites, while for U < V12, the half filled ground state gives higher preference to the state with two electrons of 
opposite spins at one site, for U > V12, the state with one electron at each site is more prefered. 

To study the transport properties through a double quantum dot system comprising of a donor and an acceptor in 
the weak coupling regime, we parameterize the different coupling strengths in the total system (system+leads) . For 
perturbation theory to be valid at temperature T, we ensure that T fc^T. More specifically in our calculations, 
we use the value of T = 0.25meV for T=300K and T = O.OlmeV for T=0.66K, which are also much smaller than the 
corresponding charging energies, e.g. Hubbard U. As our primary interest focuses on NDR effect in the system, we 
choose asymmetry in the onsite energy (Ae —e 2 -ei) to be larger than the interdot hopping parameter (t), and vary 
the Hubbard U around Ae. 

We adopt the well known exact diagonalization (ED) method to solve the Hamiltonian in equation 3 for the system 
containing two sites. As the total number of electrons, N, and z-component of the total spin, S z , commute with the 
Hamiltonian (H) and can be considered as conserved quantities, the H matrix can be diagonalized for a particular 
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FIG. 2: The number of electrons (N) in a 2 site donor-acceptor system with variation of gate bias of the system (W g ) for 
various values of the Hubbard prameter U. The inset shows the same for 17 = 5 eV with different values of nearest neighbor 
Coulomb interaction parameter, V12. Here, £2 = — £i = 2.0 eV, t = 1.0 eV and all the values of U and V12 are in eV (ft = 1). 



charge and spin sector. The Fock space then can be factored into many blocks, with largest block consisting of 4 
states with quantum numbers, number of electron (N)=2 and S z =0. The ground state energies for N=l (E\ e ), N=2 
(E 2e ) and N=3 (E 3e ) with onsite energies ei = e 2 =0 can be easily found (n = 1) 



E U = -t-Wg 

E 3e = U-t- 3W g (4) 

Thus, the gate bias window (AW S = E® e + E\ e — 2i?2 e ) over which the N=2 (half-filled) state becomes lowest energy 
state can be estimated to be, 

AW g = -2t + Va + ^{U -V 12 f + m 2 (5) 

which strongly depends on the parameters involved. E\ el E 2e and E® e are the ground state energies for le, 2e and 
3e states respectivily in the absence of gate bias. However, with inclusion of asymmetric on-site energies, the general 
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FIG. 3: Current (I) - source-drain voltage (V) characteristics of the 2 site donor-acceptor system for various values of U and 
the inset (a) represents the same for U = 5 eV with various V12 values, for £2 = — ei = 2.0eV,t — leV, T = 0.25meV and 
Temperature (T)=300K. The inset (b) shows the same, for e 2 = -ei = 2meV,t = 0.2meV, U = 4meV,T = O.OlmeV and 
Temperature (T)= 0.66K 

analytical expression for the energies become quite lengthy. For a chosen onsite energy values considering donor and 
acceptor sites, ti = — ei = 2.0eV, we plot in figure[5]the number of electrons in the lowest energy state as a function of 
gate bias. This is obtained by calculating the many-body states with minimum energy at every value of gate bias W g 
as Min [E a ). For U < AeV, with increase in onsite electron-electron interaction, there is a reduction in gate bias over 
which two electron state is stable, while on the contrary, for U > AeV, the gate bias window increases with increase in 
onsite electron-electron interaction. Furthermore, for V12 values closer or greater than U /2, the bias range over which 
two electron state is the lowest energy state, increases with increase in Vx%. This happens, because an increases in U 
by leV causes no change to E° e , an increase of leV to E® e , but an increases of more than 0.5eV to E^ e for U < 4eV 
and an increase of less than 0.5eV to E\ e for U > AeV. This is due to the fact that, for U < -ieV, two electron ground 
state gives higher preference to the state with two electrons of opposite spins at the site with lower on-site energy. 
However, for U > AeV, the it prefers the state with one electron each at the donor and at the acceptor. Hence for 
U < 4eV, an increase in U value by leV, causes an increase in the value of 2i?2 e by more than leV and an increase 
of leV to E® e , so effectively reducing the value of AW 9 . However for U > 4eV, the increase in the value of 2£ , 2 e is 
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FIG. 4: A schematic describing the transitions between the states of the donor-acceptor system, (i) in the small U regimes (ii) 
in large U regimes. The arrow direction indicates the states to which the transition occurs. 

always less than leV and hence AWj increases with increase in U values. 

For obtaining current, for every value of U and V12, the Fermi energy (Ep) is chosen as the value of the gate bias 
which ensures that two electron state is the ground state. The Fermi energy is also placed in such a way that we 
observe transition from the ground state to the state with one less electron. After fixing the Fermi energy, we have 
studied the current as a response of source-drain bias (V) in all our calculations. In figure [3j we have plotted the I-V 
characteristics of the system for a range of U and V12 values at room temperature. As can be seen clearly, low values 
of U results in step-like features in FV characterestics, while with increase in U, a rise and fall in current (a NDR 
feature) is observed for positive values of source-drain bias. Interestingly, with inclusion of nearest-neighbor Coulomb 
interaction, V12, the I-V characteristics show wide plateau region before showing NDR feature. However, the height of 
the NDR peak decreases in the positive source-drain bias region with increase in V12 . To compare our results with the 
experimental findings at low temperature, in the inset(b) of figure EH we have plotted the low temperature behaviour 




FIG. 5: The variation of (a) current (I) (b) occupation probability (Pt) of the 2e triplet state with source-drain bias, corre- 
sponding to U = 5eV and varying V12 values: Vi2=0 eV (solid line), Vi2=l eV (circle) and Vi2=3 eV (triangle). 



of I-V characteristics. Note that, the NDR peak together with the overall I-V feature 



compare fairely well with the 



17| . We also note that, there 



experimental results obtained by Tarucha et al on GaAs-based double quantum dots 
is in fact no qualitative change in the I-V characteristics except for a constant shift in bias, if we change average dot 
charge(n = 1) from one to zero in the Hamiltonian in equation 3. 



29 



30l | . However, to 



The step like feature in I-V is well understood in literature, as due to Coulomb repulsions 
understand the NDR feature, we analyze the probabiltics of occurance of various many body states. We find that 
NDR occurs when the source-drain bias drives the system from the le doublet to a higher excitation of the 2e state, 
namely the triplet states, instead of to the state with zero electron. It is because, when U is small, the ground state 
gives higher preference to the state with two electrons of opposite spins at the site with lower on-site energy. This 
allows annihilation of an electron by the electrode followed by one more annihilation leading to a transition from the 
2e singlet to the le doublet and then to the state with zero electron. However, when U increases, the ground state 
gives more weightage to the state with one electron at the donor and one at the acceptor. This allows for one electron 
annihilation from the ground state to the le doublet state, followed by a creation of an electron from the same electrode 
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to the 2e triplet state, which has the same energy as the zero electron state. Since the current at any electrode is 
calculated at steady state as the difference between the outgoing and incoming current, this transition results in a 
reduction in current leading to the negative differential resistance (NDR) peak in large U limit. A schematic figure 
describing the states involved with increase in positive source-drain bias for small and large U limits are shown in 
figure 4. Note that, with inclusion of V12, particularly for large Vn values, the charge density modulated state gets 
prominance, similar to the ground state electronic configuration as in the small U limit. Thus, with increase in V12, 
the NDR feature gets suppressed. Also, since with inclusion of V12, the gate bias range over which the 2e state 
remains the ground state differs, we pin the electrode's Fermi energy in such a way that the transition from 2e singlet 
to le doublet state occurs at the same values of V (see the inset (a) of figure 3) for a range of V12 values. However, 
with increase in positive bias, the electrochemical potential at the left electrode (hl) moves down and that at right 
electrode (/Ur) moves up, causing the transport channel e = Ex — E\ e is in resonance with the levels of the electrodes, 
where Et and E\ e are the energy levels associated with 2e triplet state and le doublet state respectively. With large 
V12 values, this channel width causes the plateau in I-V characterestics to be wider before showing NDR. 

To understand the NDR feature more clearly, and to estimate the height of the peak value in the I-V plot in figure 
3, we calculate the probability of occurance of the 2e triplet state with increase in V for a range of Hamiltonian 
parameters. In figure we have plotted the variation of current and the occupation probability of 2e triplet state 
with source-drain bias for U — 5 eV. It is clear that, when the occupation probability of the 2e triplet state starts 
increasing appreciably, the current decreases in magnitude, however, only to a nonzero value. The main point is that 
the triplet state being the blocking state suppresses the current, however, since its probability of occurance does not 
increase more than 30%, there is still some finite current (leakage current) which flows through the system. Note 
that, the I-V characteristics are asymmetric because of the inherent asymmetry in the system comprising of a donor 
and an acceptor with different site energy (e) values. 

In conclusion, we have studied transport behavior of donor-acceptor system in the Coulomb blockade regime through 
the rate equation approach. Our study shows how the variation in the on-site Coulomb repulsions can influence the 
system's response to an external source-drain bias. A strong Coulombic repulsion even results in NDR for positive 
values of source-drain bias in the I-V characteristics. Also, a strong nearest-neighbour Coulomb interaction suppresses 
the NDR like feature, taking back the system to normal Coulomb staircase regime. 
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